arXiv:1501.05498vl [astro-ph.GA] 22 Jan 2015 


Astronomy & Astrophysics manuscript no. 23 38 3 j a © ESO 2016 

August 11, 2016 


Bar pattern speeds in CALIFA galaxies 

I. Fast bars across the Hubble sequence 

J. A. L. Aguerri 1 ' 2 , J. Mendez-Abreu 12 , J. Falcon-B arioso 1,2 , A. Amorin 3 , J. Barrera-Ballesteros 1 ’ 2 , R. Cid Fernandes 3 , 
R. Garcla-Benito 4 , B. Garcla-Lorenzo 1,2 , R. M. Gonzalez Delgado 4 , B. Flusemann 5 , V. Kalinova 6 , M. Lyubenova 6 , R. 
A. Marino 7 ,1. Marquez 4 , D. Mast 8 , E. Perez 4 , S. F. Sanchez 4,9 , G. van de Ven 6 , C. J. Walcher 5 , N. Backsmann 5 , C. 
Cortijo-Fenero 4 , J. Bland-Hawthorn 10 , A. del Olmo 4 , J. Iglesias-Paramo 4,9 ,1. Perez 1112 , P Sanchez-Blazquez 13 , L. 

Wisotzki 5 , and B. Ziegler 14 

1 Instituto de Astrofisica de Canarias. C/ Via Lactea s/n, 38200 La Laguna, Spain 

2 Departamento de Astrofsica, Universidad de La Laguna, 38205 La Laguna, Tenerife, Spain 

3 Departamento de Fisica, Universidade Federal de Santa Catarina. P.O. Box 476, 88040-900. Florianopolis, S C, Brazil 

4 Instituto de Astrofsica de Andaluca - CSIC, 18008, Granada, Spain 

5 Leibniz-Institut fr Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany 

6 Max-Planck-Institut fr Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 

7 CEI Campus Moncloa, UCM-UPM, Departamento de Astrofsica y CC. de la Atmsfera, Facultad de CC. Fsicas, Universidad 
Complutense de Madrid, Avda. Complutense s/n, 28040 Madrid, Spain 

8 Instituto de Cosmologia, Relatividade e Astrofisica ICRA, Centro Brasileiro de Pesquisas Ffsicas, Rua Dr. Xavier Sigaud 150, 

CEP 22290-180, Rio de Janeiro, RJ, Brazil 

9 Centro Astronomico Hispano Aleman, Calar Alto (CSIC-MPG), C/ Jesus Durban Remon 2-2. E-04004 Almerfa, Spain 

10 Sydney Institute for Astronomy. School of Physics A28. University of Sydney. NSW 2006, Australia 

11 Dpto. de Fsica Terica y del Cosmos, University of Granada, Facultad de Ciencias (Edificio Mecenas), 18071, Granada, Spain 

12 Instituto Universitario Carlos I de Fsica Terica y Computacional, Facultad de Ciencias, 18071, Granada, Spain 

13 Departamento de Fisica Teorica, Universidad Autnoma de Madrid, Cantoblanco, E-28049 Madrid, Spain 

14 University of Vienna, Trkenschanzstrasse 17, 1180, Vienna, Austria 

Received ; accepted 


ABSTRACT 

Context. The bar pattern speed (fl b ) is defined as the rotational frequency of the bar, and it determines the bar dynamics. Several 
methods have been proposed for measuring D b . The non-parametric method proposed by Tremaine & Weinberg (1984; TW) and 
based on stellar kinematics is the most accurate. This method has been applied so far to 17 galaxies, most of them SB0 and SBa types. 
Aims. We have applied the TW method to a new sample of 15 strong and bright barred galaxies, spanning a wide range of morpho¬ 
logical types from SB0 to SBbc. Combining our analysis with previous studies, we investigate 32 barred galaxies with their pattern 
speed measured by the TW method. The resulting total sample of barred galaxies allows us to study the dependence of D b on galaxy 
properties, such as the Hubble type. 

Methods. We measured £J b using the TW method on the stellar velocity maps provided by the integral-field spectroscopy data from 
the CALIFA survey. Integral-field data solve the problems that long-slit data present when applying the TW method, resulting in the 
determination of more accurate fl b . In addition, we have also derived the ratio ft of the corotation radius to the bar length of the 
galaxies. According to this parameter, bars can be classified as fast (ft < 1.4) and slow (ft^XA). 

Results. For all the galaxies, ft is compatible within the errors with fast bars. We cannot rule out (at 95% level) the fast bar solution 
for any galaxy. We have not observed any significant trend between ft and the galaxy morphological type. 

Conclusions. Our results indicate that independent of the Hubble type, bars have been formed and then evolve as fast rotators. This 
observational result will constrain the scenarios of formation and evolution of bars proposed by numerical simulations. 

Key words. Galaxies: kinematics and dynamics - Galaxies: structure - Galaxies: photometry - Galaxies: evolution - Galaxies: 
formation 


1. Introduction 

Bars are ellipsoidal-like structures present in the inner region 
of most disc galaxies. The galaxy classification into barred and 
non-barred has been made since the pioneering Hubble morpho¬ 
logical classification. Indeed, barred galaxies constitute one of 
the branches of the so-called Hubble tuning-fork diagram (see 
iHnhhle . IT936b . 

About 40-50% of nearby disc galaxies observed in the 
optical wavelengths have been classified as barred systems 


(see, e.g., Marino va & Jo gee L 120071: iBarazza et al. L |2008; 


Agu erri et al. L l2009h . This fraction is even higher (about 60 
70%) when gala xies are observed th r ough near-infra r ed fil 


ters (see, e.g., Knapen et al. , 200(1 lEskridge et al. L [2000; 


lMenendez-E)elmestre et al. 112007 ). The bar fraction depends on 
different galaxy properties, such as Hubble type, mass, and envi- 


of the bar fraction on the Hubble type (e.r 
Buta et aL~l 120101: iBarwav et al. L [21 


, Asuerri et al. L 

2009 

Masters et al. , 

2011 


Mari nova et al. L 20121) . Several studies have recently demon- 
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strated that total stellar mass is a n important physical param¬ 
eter that regulate s bar formation (iM endez- Abreu et ah! [2010; 
iNair & Abraham L 1201 Ot ISheth et al. [ 2008 ). However, the en¬ 
vironment also plays a non-negligible role in their formation 
and evolution, as shown by Mendez-Abreu et al. (2012) when 
comparing the bar fraction in three different environments (field, 
Virgo, and Coma clusters). They show that the fraction of bars in 
the Virgo and Coma clusters are statistically different from the 
field one. They argue that bright disc galaxies are stable enough 
against interactions, whereas for fainter galaxies, interactions in 
cluster environment become strong enough to heat up the discs, 
inhibiting bar formation and even destroying the discs (see also 
ISanchez-Janssen et al. L120101) . 

Bars are very prominent features in light. In contrast, they 
represent a small fraction of the mass of the discs of galaxies. 
Hydrodynamical simulations of barred galaxies show that the 
bar component contributes on l y about 10% to 20% of the disc 
potential (see England et al. L 199 0b[ lLaine & Hell er! fl 99S 


Lindb lad et al. 11996al ; 1 Weiner et al. L 20011; lAguerri et 
20011) . Nevertheless, they produce important changes in the 


dynamics of the discs when they are present. In fact, simu¬ 
lations also show that bars are very efficient at redistributing 
angular momentum, energy, and mass on both luminous 
and dark matter galac tic comp o nents (see Weinberg! 1985 


and dark matter gaiac tic comp o nents (see Weinber g L 1 Vs 
DebatfotejfcSiellwood,, _ 1998t [20001 |Atiianaasou!a_e£ai 


2003; iM artinez-Valnuesta et al. L 20061: ISellwoodL ~ 2006; 


Sellwood & Debattista . 20061 : Villa-Vargas et al. L 20091) . 


Bars are fully characterized by three parameters: length, 
strength, and pattern speed. Both length and strength can be de¬ 
termined using optical and/or near-infrared images. In contrast, 
the bar pattern speed is a dynamical parameter, and its determi¬ 
nation requires kinematics. 

The length of the bars has be e n determined by optical vi¬ 
sual inspection iKormendv I dl979l) : lMartinl ( 1 995b. loca t ing th e 


maximum of the iso p hotal ellipticitv IWozn iak et al. l (fl 995 ); 


Marq uez et al. I ( ~999t): Laine et ahl (12002 k lMarinova &J ogee 


(2007); Aguerri et al. l (2009), or by struct ural decompositions 


2001); Aguerri et al. (2001, 2003 

20051); [Laurikainen et al. 

(20051 

2007) 

Gadotti ( 

2008); 

Laurikainen et al. (2009); 

Weinzirl et al. 

( 

2009); Gadotti 

2011). All these studies 


have shown that a typical bar radius is 3-4 kpc (e.g., 


___ typ ; _,,_. 

iMarinova & Jogeel 120071 : lAguerri et al. U2009h . There is a clear 
dependence between the bar length and other galaxy parameters, 
such as disc scale le ngth, galaxy si z e, galaxy colour, or promi¬ 
nence of the bulge lAgu erri et al. I d2005l): IMarinova & Jogeel 
d2007l) : iGadotti I d201 ll) : iHovle et al. l d201 ll) . In addition, there 
is a dependency of the bar length on the Hubble type. In 
particular, SO galaxies show la r ger bars than late-type ones 
flElmegreen & Elmegreen L 1985!: lAguerri et ah! 120091: lErwin L 
12005b Menendez-Delmestre et al. L 20071 but see also Masters et 


al. 2011). 

The bar strength is a parameter that measures the 
non-axisymmetric forces produced by the bar potential 
in th e disc of galaxies (see, e.g., lLaurikainen & Saiol 
1 20021) . The most popular methods for determining the 
bar strength are to measure the torques of the bar (e.g. . 


Combes & Sanders , 1981; Oui 

len et al. , 

2001; Laurikainen et al. , 2007; 

Salo et al. 


ifl 

10 ) 


199'. 

,201 


Buta & Block. 


, bar ellipticity 


Martinet &Friedli _ _ _ _ _ , 

d20Q2l) : Marinova&Jogee I (1 20071) : lAguerri etak _ ( 20091). 

or " 


Fouri er de composition of th e galaxy light lohta et al. I 


(1990): Marquez et al. 


lAthanassoula & Misiriotis [ 


Aguerri et al. 


[); lLaurikainen et al. 


(200C 

(20051) 


In general, SO galaxies show weaker b ars than late-ty p e ones 
(see, e.g., lLaurikainen et al. [ 120071 : lAguerri et al. L 120091 : 
iButa et alT!l2010l). 

The bar pattern speed (T2 h ) is a pure dynamical parameter 
that fully determines the dynamics of a barred galaxy. It is de¬ 
fined as the rotational frequency of the bar. Theories of stellar 
orbits have shown that Q h has a physical upper limit. Thus, a 
bar cannot extend beyond the corotation resonance (CR) radius 
(Rcr) of the galaxy. This corotation radius is the region of the 
galaxy where the angular speed of the stars of the disc in circu¬ 
lar motions equals the bar pattern speed. This limit is imposed 
by the stability of the main family of orbits forming the bar (the 
so-called x\ family; Contopoulos 1980) which are only stable 
within Rcr- Orbits in the oute r region of the disc (R > Rcr) can¬ 
not su pport a bar structure (Contopoulos , 1980; Athanas soula L 
1992). 


Measuring the bar pattern speed is technically the most 
difficult amongst the observational parameters characterising 
bars, and several methods have been developed in the literature 
to determine it. Hydrodynamical simulations of individual 
galaxies have been extensively used to determine Ob in barred 
galaxies. These methods start by computing the gravitational 
potential of the galaxy that is used to produce hydrodynamical 
galaxy models, with Qb as one of the free parameters of the 
models. Determination of the best value of Qb is done by 


and/or gas velocity field (see, e. 

g., Sanders & Tubbs , 1980; 

Hunter et al. , 1988; England et al. 

, 1990s Garcia-Burillo et al. 

19931 Sempere et al. , 1995al Lindblad et al. 

1996a; 

Lindblad & Kristen 

1996b Laine & Heller 


1999; 

Weiner et al. , 2001 

Aguerri et al. , 2001; Perez et al. , 

2004; 

Rautiainen et al. , 2008; Treuthardt et al. , 2008). Other methods 


determine the pattern speed by identifying some morphological 
galaxy features with Lindblad resonan ces. We can mention the 


positi o n of galaxy rings ( see, e . g., ButaJ, 1986; Butaetak 


see, e. g., __ _ __ _ 

1995| |Ve ga Beltran et al. L 119971: iMufioz-Tunon et al. L 2004 


spiral arms near Rcr (Canzian , 1993c Canzian & Allen , 1997; 

Puerari & Dottori, 1997; Aguerri et al. , 1998; 

Buta & Zhang 

2009); detecting the offset and shape of dust 

anes (see, e.g., 
1992); locating 

van Albada & Sanders , 1982:; Athanassoula , 


colour and star formation changes outside the bar region 
(ICeoa & Beckman! 1 1990UAguerri et al71 l2000h : or studying the 
morphology of the r esidual gas velocity field after the rotation 
veloc ity subtraction dSempere et al. L1 1995 hi iFont et alT! 1201 ll 
120141) . Nevertheless, the most accurate method for measuring 
the bar pattern speed is the non-parametric method proposed by 
Tremaine & Weinberg (1984; hereafter TW). The goal of this 
paper is to apply the TW method to a large sample of barred 
galaxies. 

The organization of the paper is as follows. The theoretical 
basis of the TW method is shown in Sect. 2. Section 3 shows the 
description of the galaxy sample. Section 4 presents the photo¬ 
metric parameters of the bars. The stellar velocity maps of the 
galaxies are shown in Sect. 5, and results are given in Sect. 6. 
The discussion and conclusions are shown in Sects. 7 and 8, re¬ 
spectively. Throughout this paper we have used the cosmology 
Ho = 70 km Mpc~* s _1 , Q,„ = 0.3, and Q,i = 0.7. 


2. Theoretical basis of the TW method 

We let (X, Y ) be a Cartesian coordinate system in the sky plane 
and take the origin at the centre of the galaxy, with X-axis coin¬ 
cident with the line of nodes (LON), which is defined as the in- 
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tersection between the sky plane and the plane of the disc of the 
galaxy. Assuming that the disc of the galaxy has a well-defined 
pattern speed G h and that the surface brightness of the measured 
tracers follow the continuity equation, the TW method is based 
on the equation: 

„ .. r_:h(Y ) r_:*(x,Y ) v WS (x,Y)dxdY {v) 

Qb sin; = - — - — - = 7^7, ( ! ) 

L 0 h(Y)f m XZ(X,Y)dXdY (X) 

where i is the galaxy inclination, Vlos (X, Y) is the line of 
sight (LOS) velocity, KA. Y) represents the surface brightness 
of the galaxy, and h(Y) is a weight function. The numerator 
of the previous equation represents the mean velocity weighted 
by the light of the galaxy ((V)). The denominator shows the 
weighted mean position of the tracers ((X)). We denote the nu¬ 
merator and denominator of Eq. (1) as ’’kinematic” and ’’photo¬ 
metric” integrals, respectively. All previous parameters can be 
measured from photometric or kinematic observations. Thus, 
i can be obtained from the ellipticity of the outermost galaxy 
isophotes, X(A, Y) can be determined from optical photometry, 
and the velocity field of the galaxy (VTos(A, Y)) can be mea¬ 
sured from long-slit or integral-field spectroscopy. In practice, 
the integrals from Eq. (1) are calculated along several directions 
oriented parallel to the LON but offset by a distance To. The 
slope of the straight line defined by the (V) versus^ (X) computed 
points represents Ob sin i (see lTremaine & WeinbergLll984l) . In 
the case of long-slit spectroscopy, the weight function is given 
by h(Y) = <5(Y - To). Similar weight functions can be defined for 
integral-field spectroscopy. Formally, integrals in Eq. 1 are over 
-oo < X, Y < +oo. Nevertheless, they can be limited to a finite 
value of X and T if the axisymmetric part of the disc is reached 
(see Sect. 6.1). 


Stellar kinematics has been used to measure Ob using the 
TW method for a large sample of early-type galaxies (e.g. . 


Kent L 19871: Merr 

field & Kuiiken , 1995; Gerssen et al. 1999; 

Debattista et al. , 

2002; Aguerri et al. , 2003F" Corsini et al. 

2003; Debattista & Williams , 2004; Corsini et al. , 2007). In 


contrast, few pattern speeds have been determined through 
this method for late-tv pe galaxies (e.g., iGerssen et al. L [2003; 
iTreuthardt et al. L 120071) . See also Corsini (2011) for a review of 
the bar pattern speed measurements using the TW method. The 
lack of measurements of bar pattern speed in late-type galaxies 
by the TW method is basically due to the effects produced on 
the determination of Ob by the uncertainties in the position an¬ 
gle (hereafter PA), i, and/or the effect of the presence of dust and 
star formation in late-type discs. Debattista (2003) demonstrates 
that errors of a few degrees in determining the PA of galaxies 
can significantly change the measured value of Ob using the 
TW method. For this reason galaxies with small errors in PA 
and/or i should be selected. In addition, the measured value of 
Ob can also be affected by dust and/or star formation in late- 
type galaxies, since in these cases the light does not trace the 
mass distribution of these galaxies. Nevertheless, these effects 
can be mitigated by computing the kinematic and photometric 
integrals of Eq. (1) using the mass distribution as weight (see 
iGerssen & Debattista L 2007 ). 

In recent years, an extension of the TW method explained 
above has been a pplied to galaxies th a t show multip l e pat¬ 
tern speeds (see i Macieiew skTl 120061 ICorsini et al. I l2007t 
iMeidt et al. L 120091) . T he TW method was al s o applied to gas 


tracers, such as CO dZimmer et al. , 2004; Rand & Wallin 

2004 

) or Her (Hernandez et al. , 2005; Emsellem et al. 

2006 

Fathi et al. , 

2007; Chemin & Hernandez , 2009 

Gabbasov et al. , 2009; 

Fathi et al., 2009). 


2.1. Fast and slow bars 

Barred galaxies are commonly classified according to the 
distance-independent ratio R = /(cr/oh, where Rqr and a h are 
the corotation and bar radius, respectively. Theoretical works 
based on stellar or bits in barred potential predict that R = 
1.2 + 0.2 (see lAthanassoula L Il992l) . This parametrization per¬ 
mits a classification of bars into fast (1.0 < R < 1.4) and slow 
ones (R > 1.4). Most observed bars have turned out to be fast 
bars; nevertheless, there a re few bars in the literature compatible 
with being slow bars (see Rautiai nen et al. L [20081) . 

There are some hints in previous works that late-type galax¬ 
ies could have higher mean values of R than early-type ones 
(see lAguerri et al.[ll998l) . More recent results based on hydro- 
dynamical models and comparing the morphology of real galax¬ 
ies with models show that while early-type galaxies always have 
fast bars, late-type gala xies host both slow and fast bars (see 
iRautiainen et al. 1120081) . Nevertheless, this result has not been 
confirmed in a significant sample of galaxies of different mor¬ 
phological types by using the TW method. The sample of late- 
type galaxies analyzed so far by the TW is very small in order 
to infer some dependence between R and the Hubble type. The 
aims here are to determine the pattern speed of a large sample 
of galaxies throughout the Hubble sequence by using the TW 
method and to study the dependence of the pattern speed of the 
bar with the Hubble type. 


3. The CALIFA sample of barred galaxies 


The barred galaxies used in the present study were taken from 
the Calar Alto Legacy Integral Field Area (CALIFA) Survey 
(ISanchez et al. L i2Q12h . CALIFA’s mother sample is formed by 
939 g alaxies selected from t he SDSS-DR7 photometric cata¬ 
logue (Abazaii an et al. L [20091) . The main selection criteria were 
angular isophotal diameter 45 < D 25 < 80 arcsec and redshift 
range 0.005 < z < 0.03. These criteria ensure that the selected 
objects fit well into the field of view (FOV) of the instrument. 
This survey aims to obtain spatially resolved spectroscopic in¬ 
formation for a fraction of the mother-galaxy sample (w 600 
galaxies), limited by available telescope time. For more details 
about the properties of the galaxies of the CALIFA mother sam¬ 
ple, see Walcher et al. (2014). This project has been the biggest 
effort with integral field spectroscopy so far. 

The CALIFA observations were carried out at the 3.5m tele¬ 
scope of the Calar Alto observatory with the Potsdam Multi 
Aperture Spectrograph (PMAS; Rot h et akl l2005h in PPAK 
mode. This instrumental mode consists of 382 fibres of 2.7 arc- 
sec diameter each. These fibres cover a FOV of 74" x 64". 
A dithering scheme of three pointings were adopted to cover 
the full FO V. This allows us to have a final resolution of 1 
arcsec (see ISanchez et al. L l2007t Perez-Galle go et al. L 1201 Ol: 


iRosales-Ortega et al. boiot Sanchez et al. [ 12012 ). Both the 
FOV and the spatial resolution of the observations make this 
dataset ideal for studies of extended disc galaxies proposed here. 

The objects were observed using two different setups. First, 
the grating called V500 shows a nominal resolution of A/AA = 
850 at 5000 A and covers from 3745 to 7300 A. The second setup 
was done using the VI200 grating with better spectral resolution 


A/AA = 1650 at 4500 A. This grating covers from 3400 to 4750 
A. In the present work we used the observations throughout the 
V1200 grating. It was selected in order to have reliable velocity 
dispersion curves of the galaxies needed for computing Rcr (see 
Sect. 5.2). The galaxy sample of the present paper was selected 
among the 200 that were first observed and reduced objects with 


3 










































































































J. A. L. Aguerri et al.: Pattern speed in CALIFA barred galaxies 


the V1200 grating. For more details about the observation strat¬ 
egy and the data reduction process, see Sanchez et al. (2012) and 
Husemann et al. (2013). 

The CALIFA mother-sample photometric properties were 
obtained from SDSS-DR7. In addition, the galaxies of the 
CALIFA mother sample were morphologically classified by us¬ 
ing the i- band SPSS ima ges by a group of five people from the 
CALIFA team (see I Walcher et al.T 120 1 4l) . The total number of 
galaxies showing strong bars (SB) on the CALIFA sample turned 
out to be 156. Only 41 of them have already been observed with 
the V1200 grating. Nevertheless, the TW method cannot be ap¬ 
plied to all galaxies. Galaxy inclination, position angle, and bar 
orientation should be considered for its application. In particu¬ 
lar, low and high inclined galaxies and those with bars near the 
galaxy’s major or minor axis should be not taken into account 
because the TW is not applicable to these objects. Therefore, we 
excluded from our analysis all bars that show PA less than 10° 
with respect to the major or minor axis of the galaxies and either 
those face-on (b/a > 0.86) or edge-on (b/a < 0.34) galaxies. 
This restriction reduces the total number of barred galaxies to 
20. In addition, five more galaxies were excluded for other rea¬ 
sons, such as the large number of field stars, the non-flat rotation 
curve at large radius, or/and the insufficient quality of the stellar 
velocity maps. 

The final selected sample of barred galaxies contains objects 
that cover the morphological types: 1 SB0, 2 SBO/a, 1 SBa, 1 
SBab, 6 SBb, and 4 SBbc. Similar to the CALIFA mother sam¬ 
ple, our galaxy sample is dominated by spirals (SBb-SBbc). This 
is appropriate for filling the gap of pattern speeds measured by 
the TW method in late-type galaxies. 

Figure Q] shows the comparison of the distribution of mor¬ 
phological types, redshift, absolute r-band magnitudes, and bar 
radius between the barred galaxies in the CALIFA mother sam¬ 
ple and the selected barred galaxies of the present work. Our 
sample is a small subsample of the total bar population of 
the CALIFA mother sample. Nevertheless, the Kolmogorov- 
Smirnov test indicates that our subsample does not have a sta¬ 
tistically different distribution of morphological types, redshift, 
and absolute magnitudes from the total sample of barred galax¬ 
ies. The actual sample loss bars have radii larger than 10 kpc. 
The lack of large bars is due to: i) galaxy inclination, ii) bar ori¬ 
entation, and/or iii) no observations using the V1200 grating. We 
never cut the bar sample according to bar size. It is also worth 
noticing that the galaxies presented here are all bright galaxies 
( M r > -19.0). In this range of magnitudes, CALIFA is represen¬ 
tative of the local Universe. No dwarf galaxies have been con¬ 
sidered in the present work. The r-band SDSS images of the 
final sample are shown in Fig. [2] and the main parameters of the 
galaxies are given in Table [T] 

In summary, the sample presented in this work overcomes 
two major problems of other samples presented in the literature. 
First, integral-field observations such as those presented here get 
rid of many of the problems of long-slit observations. Second, 
our sample is dominated by late-type galaxies, the missing piece 
of information in the TW studies of bar pattern speeds. 

4. Photometric parameters of the barred galaxies 

4.1. Inclinations and position angles 

Analysis of the galaxy isophotes provides the inclination^, 
PA, and other important information about the different struc- 

1 The inclination ( i ) of a disc is given by b/a = (q 1 2 + (1 - q) 2 x 

cos 2 (i)) 111 , where b/a is the ratio between the observed minor and major 






Fig. 1 . Normalized distribution of morphological types (a), red¬ 
shift (b), absolute r-band magnitudes (c), and bar radius (d) of 
the CALIFA’s mother-sample barred galaxies (full black line his¬ 
tograms) and the barred galaxies selected in the present paper 
(dashed blue-line histograms). 


tural compone n ts of g alaxies (see, e.g., IWozniak et al. L 119951: 
Ague rri et al. L l2000 bl. The galaxy inclination makes the 
isophotes of the galaxies appear as ellipses in the plane of the 
sky. We fitted the isophotes of our galaxies by ellipses using the 
ELLIPSE routine from the IRAF package dJedrzeiewski L ll987t. 
The ellipticity and PA isophotal radial profiles obtained from 
these fits are shown in Fig. [3] 

The ellipticity radial profile of an unbarred spiral galaxy 
grows from almost zero values at the centre of the object up to 
a constant value at large radii, which corresponds to the galaxy 
inclination. The PA isophotal radial profile also reaches a con¬ 
stant value at large radii, which corresponds to the orientation 
of the LON of the galaxy. This behaviour can be observed in 
the ellipticity and PA isophotal radial profiles shown in Fig. [3] 
We therefore computed the i and PA of our barred galaxies by 
averaging their outer isophotes. The results are given in Table Q] 


4.2. Determination of the bar radius 


The determination of the bar radius is not an easy task. Several 
methods have been proposed in the literature during the past few 
decades. We used three of the most popular methods (ab,i.flb, 2 , 
and flb, 3 ) for estimating the bar radius (a h ) of our galaxies. 

Two of the bar radius measurements (ab.i and flb, 2 ) were ob¬ 
tained by using the information of the ellipticity and PA isopho 
tal radial pr ofiles and the di f ferent methods proposed in the lit 
erature fsee iMarauez et al. L 1999t lAthanassoula & Misiriotis 
120021: iMichel-Dansac & WozniakTl2006l: lAguerri et al. L l20(f ~ 


These methods are based on the peculiar features produced by 
the shape and orientation of the stellar orbit s of barred galax¬ 
ies dContopoulos , 1980: lAthanassoulall 19921) . In particular, the 
galaxy isophotes are almost circular at the galaxy centre, because 
of either seeing effects or the presence of a spherical bulge. As 
we get away from the centre, there is a general increase in the el¬ 
lipticity up to a local maximum, and then it suddenly decreases 
towards a minimum at the location where the isophotes become 


axis of the disc, and q represents its intrinsic thickness. In this paper 
we have assumed the thin disc approximation (q = 0). Thus, we have 
computed the inclination of the galaxies by cos(i ) = b/a 
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Table 1 . Main parameters of the CALIFA barred galaxy sample 


Galaxy 

RA (J2000) 
(hh:mm:ss) 

DEC (J2000) 

(° 

Morph. Type 

i 

(degrees) 

PA 

(degrees) 

Re 

(arcsec) 

m r 

(mag) 

NGC 0036 

00:11:22.3 

06:23:22 

SBb 

57.2±3.7 

23.4 ±1.3 

14.9 

12.7 

NGC 1645 

04:44:06.4 

-05:27:56 

SBOa 

64.5±0.6 

84.7 ±0.4 

9.2 

12.7 

NGC 3300 

10:36:38.4 

14:10:16 

SBOa 

57.2±0.7 

172.0±0.3 

12.0 

12.3 

NGC 5205 

13:30:03.6 

62:30:42 

SBbc 

50.0±0.7 

170.1 ± 1.6 

15.9 

12.5 

NGC 5378 

13:56:51.0 

37:47:50 

SBb 

37.8±3.1 

86.5 ±5.4 

19.9 

12.5 

NGC 5406 

14:00:20.1 

38:54:56 

SBb 

44.9±0.4 

111 .8±0.8 

18.1 

12.2 

NGC 5947 

15:30:36.6 

42:43:02 

SBbc 

44.6±1.4 

72.5 ±3.1 

12.4 

13.4 

NGC 6497 

17:51:18.0 

59:28:15 

SBab 

60.9±0.7 

112.0±0.4 

10.7 

13.0 

NGC 6941 

20:36:23.5 

-04:37:07 

SBb 

42.3±0.8 

127.5±2.7 

14.3 

13.1 

NGC 6945 

20:39:00.6 

-04:58:21 

SB0 

51.3±2.3 

126.1 ± 1.3 

9.3 

12.4 

NGC 7321 

22:36:28.0 

21:37:19 

SBbc 

48.3±0.5 

13.4 ±1.3 

12.0 

12.9 

NGC 7563 

23:15:55.9 

13:11:46 

SBa 

55.8±2.4 

149.8± 1.3 

8.1 

12.4 

NGC 7591 

23:18:16.3 

06:35:09 

SBbc 

57.6±0.5 

144.0±2.0 

12.0 

12.6 

UGC 03253 

05:19:41.9 

84:03:09 

SBb 

56.8±1.6 

92.0 ±1.7 

11.8 

13.2 

UGC 12185 

22:47:25.0 

31:22:25 

SBb 

64.0±0.6 

161,0± 1.0 

8.8 

13.5 


- Note: Columns are: (1) NGC name of the galaxy; (2) galaxy right ascension, (3) galaxy declination, (4) morphological type, (5) galaxy 
inclination measured from isophotal ellipticity profiles on i-band SDSS images, (6) galaxy PA measured from isophotal PA radial profiles on 
i-band SDSS images, (7) effective radius in r-band from SDSS-DR9, (8) model r-band magnitude from SDSS-DR9 


axisymmetric (disc region) in the face-on case. Typical mea¬ 
surements of the bar radius using the ellipticity profile involve 
measuring the position of this maximum and minimum. In fact, 
they represent two extreme cases (Michel-Dansac & Wozniakl 
120061) . and therefore they can be understood as the lower and 
upper limits to the bar radius measurement. We adopted the po¬ 
sition of the local maximum of the isophotal ellipticity of our 
galaxies as one measurement of the bar radius of our galaxies. 
Table [2] shows the values of the bar radius for our galaxies a\, \ 
determined by this method. 

The second estimate of the bar radius (a h 2 ) takes the infor¬ 
mation provided by the isophotal radial PA profiles into account. 
The PA radial profile is also characteristic of barred galaxies, 
since it is constant in the bar region and then changes to fit the 
outer d isc orientation (e.g., Wozn iak et al. Lfl995L I Aguerri et al. L 
2000b). A typical bar length is measured at the radius where the 
position angle changes by APA with respect to the value corre¬ 
sponding to the maximum ellipticity. Usually APA=5° is a good 
value for the bar radius (see e.g., [Aguerri et al. Il2009t) . The value 
of the bar radius obtained with this method could be correlated 
to a/, i. This is due to the behaviours inside the bar region of 
the ellipticity, and PA profiles are similar among galaxies. The 
values of « b .2 for our galaxies can be seen in Table[2] 

Fourier decomposition has been used extensively in charac¬ 
terising structures, like bars, which represents a bisymmetric de¬ 
parture from axisymmetry. We used the method based on Fourier 
decomposition of the light distribution of the galaxies proposed 
by Aguerri et al. (2000). Following this method, the bar radius is 
determined by the ratios of the intensities in the bar and inter-bar 
regions. The azimuthal surface brightness profiles of the depro- 
jected galaxies were decomposed in a Fourier series. The bar 
intensity, 4, is defined as 4 = 4 + 4 + 4 + 4 (where 4,4,4 
and 4 are the m = 0,2,4, and 6 terms of the Fourier decompo¬ 
sition, respectively). Similarly, the inter-bar intensity is defined 
as 4> = 4 - 4 + 4 - 4- The bar region is defined as the region 
where 4/fib > 0.5 x [majc(4/4>) - min{ 4/fib)] + min(I\,/I,b). 
The bar radius (a/,^) is identified as the outer radius at which 
4/fib > 0.5 x [max(4/fib) -wtfn(4/fib)] +wm(4/fib)- Numerical 
simulations have shown that this method determines the bar ra¬ 
dius within 8% accuracy except for very thin bars (see, e.g., 


lAthanassoula & Misiriotis 11200 21. The values of for our 
sample are shown in Table [2] 

Aguerri et al. (2009) demonstrated that the best method for 
determining the bar radius depends on the shape of the sur¬ 
face brightness profile of the bar (see their fig. 4). This bar 
type is only available after accurate multi-component surface- 
brightness decomposition of photometrical images. They also 
demonstrated that ab.i < c'b ,2 < «b ,3 independent of the bar type. 
In addition, the real bar radius is always between «b,i and a h . 3 - 
For these reasons, we consider a\, as the mean of flb.i,flb, 2 , and 
ab, 3 - Moreover, the upper and lower uncertainties of the bar ra¬ 
dius are the differences between a h and a h ,i, Ob, 2 , and c/ b 3 , re¬ 
spectively. 

5. Stellar velocity maps 

The stellar kinematics of the galaxies was measured from the 
spectral datacubes observed with the V1200 grating. The full 
description of the procedure will be explained in Falcon-Barroso 
et al. (in prep.). For the sake of clarity, we briefly summarise the 
process here. 

In the first step, the spaxels of the datacube were Voronoi- 
binn ed to achieve a limi t ing si gnal-to-noise ratio S/N > 20 
(see ICapnellari & ConirTl 1 20031) . while spectra with S/N < 3 
were not considered. The values of the LOS velocity and ve¬ 
locity dispersion were obtained by fitting the binned spectra us¬ 
ing the penalised pixel-fitting method (pPXF) from Cappellari & 
Emsellem (2004). The fit takes the continuum and the galaxy ab¬ 
sorption features presented in the wavelength range of the spec¬ 
tra into account. Emission lines were masked where present in 
the spectra of the galaxies. A non-negative linear combination 
of a_subset of 328 stellar templates from the Indo-US library 
(IValdes et al. L12004) were used for the fit of the spectra. This 
subset was carefully selected in the parameter space defined by 
the stellar properties: T e ff, log(g), and [Fe/H]. Errors in both ve¬ 
locities and velocity dispersion were obtained via Monte Carlo 
simulations. 

The stellar velocity maps of the galaxies considered in the 
present paper are shown in Fig [4] In these maps, the systemic 
velocity of the galaxies was computed as the average velocity of 
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Table 2. Bar de-projected radius measurements for the CALIFA 
barred galaxy sample 


Galaxy 

a b,l 

a b,2 

a b,3 

flb 


(arcsec) 

(arcsec) 

(arcsec) 

(arcsec) 

(1) 

(2) 

(3) 

(4) 

(5) 

NGC 0036 

20.8 

24.7 

15.1 

20.2«-; 

NGC 1645 

13.6 

17.3 

17.0 

16.0 

NGC 3300 

13.4 

17.6 

17.4 


NGC 5205 

14.9 

18.3 

19.8 

17.7+18 

NGC 5378 

25.3 

34.4 

23.4 

'jn '7+4.3 

A ' 6 7 

NGC 5406 

20.0 

23.2 

20.0 

21 1 + ^ 

21 

NGC 5947 

9.6 

12.5 

10.7 

10 - 9 -' 

NGC 6497 

12.6 

15.6 

16.0 

14.7Z- 

NGC 6941 

14.6 

18.0 

15.0 

15.9*2 j 

NGC 6945 

13.8 

16.4 

19.2 

16.5 + ^ 

NGC 7321 

10.4 

14.1 

11.7 

12A H 

25 . 5/54 

NGC 7563 

22.8 

30.9 

22.7 

NGC 7591 

11.0 

15.0 

14.5 

13 5 +ls 

UGC 03253 

14.5 

18.0 

14.9 

UGC 12185 

12.0 

24.3 

24.5 

20.3/f, 


Notes. Col. (1) Galaxy name. Col. (2): bar radius measured from the 
maxima of the ellipticity radial profile. Col. (3): bar radius measured 
when the PA in the bar region changes more than 5°. Col. (4): bar ra¬ 
dius measured from the Fourier bar/interbar intensities. Col. (5): mean 
of the three bar-length estimations. Errors correspond to the maximum 
differences between the mean and the three measurements. 

the stars in the central 5" aperture and subtracted from the veloc¬ 
ity maps. The stellar velocity obtained for our galaxies was used 
for computing the kinematic integrals of the TW method (see 
Eq. 1). Figure [5] shows stellar-streaming mean velocities with 
V sys subtracted, which were obtained along the galaxy’s photo¬ 
metric major axis. 

6. Results 

6.1. The bar pattern speed 

The bar pattern speed of the galaxies was measured by applying 
the TW method as described in Eq. (1). This method requires de¬ 
termining the mean weighted position ((X)) and velocity ((V)) 
of stars in several slits oriented along and offset with respect 
to the LON of the galaxies. We refer to these slits defined in the 
integral-field datacubes as pseudo slits. Depending on the bar ra¬ 
dius and its orientation with the LON, we have used three or five 
pseudo slits of 1 arcsec width each and a minimum separation of 
2 arcsec between them to avoid repeated information. Light and 
mass were used as two different weights for computing (X) and 
( V ) for each pseudo slit. 

The high S /N ratio of the broad-band photometric images 
has led in the past to their being broadly used for computing 
the photometric integrals of the TW method. Nevertheless, the 
trace of different stellar populations between the photometric 
and spectroscopic data, caused by different wavelength coverage 
and problems with the positioning of the pseudo slits in the pho¬ 
tometric images, can affect the computation of the photometric 
integrals. These problems can be solved by using integral-field 
data. In particular, we computed the mean position of stars along 
the pseudo slits by using the surface brightness distribution ob¬ 
tained directly from the CALIFA datacubes of the galaxies. The 
surface brightness map of each galaxy was obtained by sum¬ 
ming up all the flux from each spectrum of the datacube in a 
wavelength window of 150 A width, and centred at 4575 A. 


This wavelength window was used because no prominent 
emission lines are observed in this range. It is also important to 
notice that the light distribution does not always trace the mass 
distribution of the galaxies. This is especially true in late-type 
galaxies where star formation is common. To get rid of this, we 
also computed the mass-weighted mean position of stars by us¬ 
ing the profiles extracted along the pseudo-slit positions in the 
mass maps of the galaxies provided by Gonzalez Delgado et al. 
(2013) and Perez et al. (2013). Thus, for each galaxy we have 
two values of the photometric integrals of the TW: the mass- 
and light-weighted photometric integrals. 

Two different approaches were used to measure the weighted 
LOS stellar velocity (( V )) for each pseudo-slit position. The first 
approach was by computing the integrals in the numerator of Eq. 
(1) using pseudo slits placed directly in the velocity maps of the 
galaxies. The second approach consists in a weighted sum of the 
raw spectra from the datacube along the pseudo slits. This results 
in a new, single spectrum for each pseudo slit, which was then 
analysed using the pPXF method as explained in the previous 
section. Thus, (V) is the radial velocity obtained from the fit to 
this single spectrum. In the two previous approaches, (V) was 
obtained by also using light and mass as weights. In all cases, 
Monte Carlo (MC) simulations were used to compute the errors 
of <V>. 

Formally, all the TW integrals shown in Eq. (1) are over 
—oo < X < oo. But, they can be limited to -Aj nax < X < X lmx 
if X max reaches the axisymmetric part of the disc. We used the 
maximum value of allowed by our velocity maps. In most 
of the cases, this value was X m!lx = 20-30 arcsec (see Fig. 0. 
For our galaxies we have a mean value of A r max /7zdi sc =L6, where 
/zdisc is the exponential scale length of the disc of the galaxies. 
The disc scale length was determined by fitting the outermost re¬ 
gions (outside of the bar region) of the r-band isophotal surface- 
brightness profiles of the galaxies by an exponential law (see 
iFreemanl 1 1970 1. 

The pattern speed of the galaxies was determined by the 
slope of the straight line fitted to the (V) vs (X) points. The un¬ 
certainties in Qb were obtained by using MC simulations tak¬ 
ing the uncertainties in the PA of the galaxies into account. The 
PA uncertainties were distributed according to Gaussian distri¬ 
butions. We have four different determinations of the pattern 
speed (Qb.i to Q h4 ) of the galaxies depending on both the weight 
used for computing the TW integrals and the method used for 
the kinematic integrals. The TW kinematic integrals were ob¬ 
tained by computing the integrals shown in the numerator of Eq. 
(1), using the velocities provided by the galaxy velocity maps 
along pseudo-slit positions for Ob.i (light-weighted) and Db ,3 
(mass-weighted). The sum of the spectra from the galaxy dat¬ 
acube along the pseudo-slit position was used for computing 
the kinematic integrals in Qb ,2 (light-weighted) and Q h4 (mass- 
weighted). Figure[6]shows the (V) vs (X) line fits for the sample 
galaxies. The linear fits shown in Fig. [6] do not take the uncer¬ 
tainties of the photometrical integrals into account. The errors 
of the photometric integrals are in all cases smaller than 0.1 ex¬ 
cept for one pseudo slit in the case NGC0036, which is 0.14. 
Table 0 shows the values of Q h for the CALIFA barred galaxies 
presented here. 

Figure [7] shows the comparison of the values of Q b sinz ob¬ 
tained by different methods. In general, no systematic and/or sig¬ 
nificant differences in the value of Qb sin i can be seen. Only 
three galaxies (NGC5406, NGC5947, and NGC6497) present 
some value of Q/, that does not agree, within the errors, with 
the others. It is not clear to us which is the main reason for this 
difference, since these galaxies do not present irregular features 
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Fig. 2. SDSS r-band images of the CALIFA barred galaxies presented in this study. In all images, north is up and east is left. The full 
and dashed lines represent the considered slits for measuring the kinematic and photometric integrals. The full line also represents 
the line of nodes of the galaxy. 


in their rotation velocity maps (see Fig. [4]) and/or velocity curves 
(see Fig. 0. In addition, they are galaxies of different morpho¬ 
logical types, which seems to indicate that these differences are 
not related to different amounts of dust, gas, or star formation in 
their discs. These differences could simply reflect the observa¬ 
tional uncertainties in the determination of Cl/,. 

6.2. Determination ofR- /?cR/«b 

The dimensionless quantity R = Rail nr is defined as the ratio 
between the corotation radius (Rcr) of the galaxy and the bar 
radius (an). Unlike the bar pattern speed, measure 'R requires 


some modelling to recover the rotation curve from the observed 
stellar streaming velocities. 


The corotation radius is derived from the circular velocity 
of the galaxy. To obtain the circular velocity ( V c ) from the ob¬ 
served stellar streaming velocity (V*), the asymmetric drift cor¬ 
rection is needed. One of the aims of the present work has been 
to compare our pattern speeds for SBb and SBbc galaxies with 
those from early-type systems from the literature.To keep the 
homogeneity in this comparison, we followed the same approx¬ 
imation for the asymmetric drift correction as in previous works 
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Fig. 3. Isophotal ellipticity and PA radial profiles from ellipse-fitting of the barred galaxies of the sample. The horizontal lines show 
the measured ellipticity and PA corresponding to the disc. 


Table 3. Pattern speed and corotation radius of the CALIFA barred galaxies 


Galaxy 

flb.l 

Rcr.i 

f\,2 

^CR,2 

fib, 3 

RcR,3 

Hb,4 

^CR,4 


(km s -1 arcsec -1 ) 

(arcsec) 1 

^km s -1 arcsec -1 ) 

(arcsec) 

(km s -1 arcsec -1 ) 

(arcsec) 

(km s -1 arcsec -1 ) 

(arcsec) 

NGC 0036 

13.2±3.7 

16-91“ 

17.4+5.2 

12.6 + “ 

19.3±4.3 

■W 

19.5±4.3 

TTT^ 

NGC 1645 

-10.9±11.6 

18.9i 5 , 5 

-22.4±9.4 

12-1 ft 

-13.7± 12.4 

19 l lf - 3 

Lyi 5.3 

-20.2+11.7 

13-it 

NGC 3300 

-9.8±3.1 

23.2 + ”| 

-9.0±2.4 

24.5 « 





NGC 5205 

17.1 ±4.3 

10-22.4 

15.1 ±2.8 

ii4| 

11.5±2.8 

H- 2 K 

13.8+1.7 

12-3|;I 

NGC 5378 

-11.2±6.1 

16 . 31 ‘ 3 6 ° 

-9.8±4.5 

18.0+L, 0 

-9.3±4.5 

17.51 15 7 ° 

-9.0+3.9 

1 q 7+12.4 
-7.8 

NGC 5406 

16.3±9.3 

15.51^,° 

22.8±8.0 

11 . 01 ^ 

7.2±5.6 

14- 7 8.7 

9.1+ 3.8 

97 4+14.7 
^' * -8.3 

NGC 5947 
NGC 6497 

-16.4±7.3 

-32.4±5.5 

11 6 +g3 

11 ‘ u -4.8 

7-4^ 

-31.7±4.2 

-42.7±7.4 

5 -8l!3 

5 - 5+2 ii 

-18.5±4.9 

-21.4±3.6 

10.7+11 

8-7lff 

-16.4+5.6 

-26.6+2.7 

11-2111 

8 . 811-2 

NGC 6941 

12.8±10.4 

14 ’ 6- vC 

18.4±9.5 

10 . 7 ^ 

9.4±9.3 

14 .HS 9 

21.9+6.2 

9 2 +31 
y -2.2 

NGC 6945 

-10.3±3.6 

20.21” 

-15.5±2.1 

13-0 ± 

-9.7±2.4 

20.3™ 

-13.6+1.6 

14 7+3.0 
-2.6 

NGC 7321 

-16.7±3.8 

15 6 +3 - 7 

aj - u _4.4 

-21.6±6.5 

11-8^o 

-18.7±4.0 

14.61|;| 

-16.7+4.3 

15 . 312 I 

NGC 7563 
NGC 7591 

7.6±3.3 

-9.8±4.6 

OO 0 + 18.4 

JJm - -10.2 
I8.5I52 

4.5±2.7 

-13.9±5.0 

4c 4 + 19.6 

1 9 0+6.6 
LJ ‘^-4.3 

8.0±3.0 

-10.5±3.5 

09 9+13.5 
-8.4 

18.314° 

3.9+2.1 

-7.7+4.4 

ci 0 + I6.O 

J1 * — 17.1 
20.6 + )“ 

UGC 03253 

-10.5±3.1 

18.0 + a 1 

-15.5±3.1 

11. 9^4 

-12.8±2.4 

16.21|| 

-12.6+2.3 

14.6111 

UGC 12185 

13.3±1.7 

i 6 - 7 i4 

10.0±2.0 

21.6q 1 

12.8±1.6 


11.6+1.7 

18.4 +10 fi 1 


- Note: NGC 3300 has no values of and £2^4 because the mass distribution of this galaxy is not available. 


(see Debattista et aI~ll2002H Aguerri et al, L 

120031) . Thus, the asymmetric drift equation 




becomes 


Vl = V: + 


cr: 


obs 


sin 2 /(I + 2a 2 cot 2 i) 


/ 1 

2 \ 

2R — + — - 1 

\Rd 

R<r) 


( 2 ) 


where R, V t , and cr 0 i, s are the distance to the galaxy centre, 
the LOS velocity, and velocity dispersion along the major axis 
of the galaxy, respectively. Parameter a is defined as the ratio 
between the perpendicular and radial disc velocity dispersions. 
The value of a can vary from early- to late-type galaxies. We 


have used the values given by Gerssen & Shapiro (2012) for the 
different morphological types of galaxies. In particular, we have 
used a = 0.85 + 0.15 for SBO-SBO/a, a = 0.86 ± 0.24 for SBa- 
SBab, and a = 0.62 ± 0.2 for SBb-SBbc. In addition, Rj and 
are the scale lengths of the surface brightness and velocity dis¬ 
persion pr ofiles, assuming an exponential law for the disc (e.g., 
lFreemanlll970l) . 

This correction was applied to the velocity data along the 
major axis of the galaxy discs and beyond to the bar region. 
The value of V Ci fl at for the galaxies was obtained by averag- 
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Fig. 4. Binned version of the stellar velocity maps of the barred galaxies in our sample. The orientation of the maps is as in Fig. 
[2] Blue colours represent approaching velocities, while red ones show receding velocities. The full and dashed lines represent the 
considered slits for measuring the kinematic and photometric integrals. The full line also represent the line of nodes of the galaxy. 


ing V c outside of the bar region. The Tully-Fisher relation (TF; 
ITullv & Fisher L1197 7) shows that there is a correlation between 
the circular velocity of the galaxies at large radii and their abso¬ 
lute magnitudes. We located our galaxies in the TF relation by 
using the V Ci fl at obtained from the asymmetric drift correction. 
Figure [8] shows the TF relation for our galaxies and other spi¬ 
ral galaxies obtained from the literature (see lReves et al. LfeOl If) . 
We notice that our galaxies are located within the noise of the TF 
relation. This indicates that the Vc.flat we obtained is consistent 
with what is proposed by the TF relation for spiral galaxies. 

In addition, we found in the literature the maximum cir¬ 
cular velocity measured from HI data for five of our galax¬ 


ies (see iTheureau et al. I 119981) . In all cases, this HI rota¬ 
tion velocity agrees within the errors with our values of 
Vcjiat■ For these five galaxies, the ratios between the maxi¬ 
mum circular HI rotation velocity ( V c hi) and our V c fi at are 
Vc.mlVc.fia, = 1.12,1.14,0.99,0.88, and 1.05 for NGC0036, 
NGC6941, NGC7321, NGC7563, and UGC03253. This indi¬ 
cates a maximun uncertainty of 14% in our flat rotation deter¬ 
minations. 

The value of the corotation radius is then given by Rqr = 
Fc.flat/^b- This assumes that the rotation curve of the galaxies 
is flat. The values of Rqr for our galaxies are given in Table [3] 
We finally measure R as Rck/u^. The values of R are given in 
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Fig. 5. LOS-streaming velocities of the galaxies along their major axis. 


TablelU Errors in Tables[3]and[4]were computed by Monte-Carlo 
simulations. 

6.3. Fast or stow bars 

The galaxy sample presented in this work complemented by 
those available in the literature (see Table 1 from Corsini 2011) 
make a total of 32 galaxies with bar pattern speed measured by 
the TW method. This represents an appropriate sample for infer¬ 
ring strong observational constraints on the value of R in bright 
barred galaxies. Table [4] shows that all bars have values of R, 
within the errors, compatible with being fast bars (R < 1.4). We 
computed the probability of having R > 1.4 in the whole sam¬ 
ple by using Monte Carlo simulations that take the quoted errors 


into account. We cannot rule out a fast bar in any of our galaxies 
(at 95% probability). 

We have also determined the mean value of R for the total 
sample of 32 galaxies with fib determined by the TW method. 
Taking our four different measurements of R into account, we 
obtained < R x >= 1.2!°’, < R 2 >= 1.0+°;’, < R 3 >= 1.2+°;’, 
and < R 4 >- 1.1!°’. Table 4 shows that some of our values 
of R have large errors. We selected a subsample of galaxies 
with small uncertainties in R. In particular, we considered those 
galaxies from the literature and from our measurements with un¬ 
certainties smaller than 30% in determining R. For this subsam¬ 
ple of galaxies, we have < R\ >= l.O!^, < R 2 >= 0.9!°2, 
< R 3 >= 1.0!°i and < R 4 >= O^!^. In all cases, the mean 
values indicate that bars finish near corotation in agreement with 
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Fig. 6. Values of (V) vs (X) for the galaxies of the sample. The different colours represent the different techniques used for computing 
the kinematic and photometric integrals of the TW method. Black and green lines represent light-weighted photometric integrals, 
with kinematic integrals computed using velocity maps or summing spectra along pseudo slits, respectively (corresponding to £2b,i 
and Ob. 2 ), Red and blue lines represent mass-weighted photometric integrals with kinematic integrals computed using velocity maps 
or summing spectra along pseudo slits, respectively (corresponding to Qb ,3 and T2 h 4 ). 


numerical simulations of barred galaxies (see lAthanassoula L 
1199211 . 

6.4. Dependence of the bar pattern speed on the 
morphological type of the galaxy 

Figure |9]shows the variation in R with the Hubble morphological 
type for the 32 galaxies with Ob determined by the TW method. 
We also computed the mean values of R for galaxies in three 
Hubble-type bins SBO-SBO/a, SBa-SBab, and SBb-SBbc. Table 


[5]shows these mean values. The mean values of R for SBb-SBbc 
galaxies are always lower than for early-type ones. Nevertheless, 
this trend is not significant when considering the errors. 

As previously noted, the errors showed in Tableware large. 
We investigated whether these large errors in R are masking a 
trend with Hubble type. To analyse this possibility, we have also 
studied the trend in R with the Hubble type for the subsample of 
galaxies with uncertainties smaller than 30% in the determining 
R. In this case, we have only compared SBO-SBO/a and SBb- 
SBbc Hubble types, because we have no galaxies in the morpho- 
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Fig. 7. Comparison of the bar pattern speeds of the galaxies using different methods. 
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M r (mag) 

Fig. 8. Tully-Fisher relation of our CALIFA galaxies (black sym¬ 
bols) and spiral galaxies from Reyes et al. (2011) (grey sym¬ 
bols). The circular velocity of the CALIFA galaxies was op- 
tained from the asymmetric drift correction applied in this study. 
The full line represent the fit of the TF relation for the galaxies 
from Reyes et al. (2011). 


logical bin SBa-SBab. Table [6] shows these values. In this case, 
the errors are considerably smaller. Although the mean values 
of R are also systematically smaller for late-type galaxies, these 


Table 4. Dimensionless quantity R for the CALIFA barred 
galaxies 


Galaxy 

% 

Ri 

r 3 

7?4 

NGC 0036 

0.9 +a; * 

0 .6!- 

0 6 +u ' 2 

u 0.1 


NGC 1645 
NGC 3300 

1 3 +A - 7 
isM 

1 -0.4 

0 -8!; 

1 6 +0 - 5 

1 ‘ -0.4 

11 +Y:1 

1 * 1 -0.4 

0-8 il 

NGC 5205 

o-6!:r 

0 7 +0 - 2 
-Q.l 

0 q+ 0-3 
u -0.2 

0 7 +ai 
-Q-i 

0 7 +0 - 4 

w -0.3 

NGC 5378 

0 -6!| 

0 6 +0 - 4 
u -0.2 

0 6 +0 - 5 

VJ ‘ -0.3 

NGC 5406 

0 7 +06 
u -0.2 

0 5 +a2 

w -Q.l 

l 4+1.5 

1 ' -0.6 

1 3 +av 

1 -0.4 

NGC 5947 

1 o +as 

04 

0 5 +a2 

0 9 +0 - 5 

y 04 

1 0 +0 ' 6 

04 

NGC 6497 
NGC 6941 
NGC 6945 
NGC 7321 

0 4 + °-* 
09^ 

7,5*3 
L3+i 

0 3 + ^ t 
a«) 

L0 n. 

0 6 + °^ 

1^:1 

IV® 

03 


NGC 7563 

1 - 2 -0 3 

1 9 ^ 

1.7 qj 

2 3 + ^ 

T'^3 

09^ 

"■ y -Q2 

NGC 7591 
UGC 03253 

1 4+Y.o 

i^ 

1 * -0.3 

or*2 
u - 7 -0.2 

i a+0.6 

09^ 

u, *-0.2 

UGC 12185 

0 9 +a4 

v - y -Q3 

i r )+0.6 

1 -0-5 

0 9 +0 - 4 

W, ^-0.4 

i o +0 - 5 

a - u -0.4 


differences are again compatible within the uncertainties. We can 
therefore conclude that no significant trend is observed in the 7? 
values of early and late-type galaxies. 

For the sake of comparison in Fig. [9] we have also included 
the mean values of R in three different morphological ranges 
(SBO+SBO/a, SBa+SBab, and SBb+SBbc) obtained through hy- 
drodynamical simulations by Rautiainen et al. (2008). In contrast 
to the results shown by Rautiainen et al. (2008), the application 
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Table 5. Mean values of R for different Hubble types. 


Hubble type 

Ri 

Ri 

r 3 

Ra 

Ngal 

SBO-SBO/a 

l -5+0.6 

1 -0.5 

1 9+0.6 

A - -0.4 

1 -5+0.6 

1 -0.5 

1 9+0.6 

1 ■ -0.4 

17 

SBa-SBab 

i 9+O.6 
1,z '-0.5 

l-4lVo 

1 2 +a * 

1 ‘ -0.5 

1 3 +1 ° 

1 -0.8 

3 

SBb-SBbc 

1 1+0-6 
04 

°- 8 t 

1 0 +os 
u 03 

n g+o.6 

y, - y 03 

12 

All types 

1 2 + ^ 

1 -0.5 

1 0 + ^ 

1 -0.4 

1 2 +ij:7 
-0.4 

l i +!i:7 

1 -0.4 

32 


- Note: N ga i represents the number of galaxies for each morphologi¬ 
cal type bin. 

Table 6. Mean values of R for different Hubble types for the 
subsample of galaxies with smaller uncertainties in the determi¬ 
nation of R. 


Hubble type 

Ri 

Ri 

r 3 

Ra 

Ngal 

SBO-SBO/a 

1 2 +{>:i 

-Q.2 


1 2 +u - 2 * 

-Q-2 

L1 -S4 

13 

SBb-SBbc 

°- 9+ oJ 

o-r 

0 9 +0 - 5 

yj.y q| 

o.r 

10 

All types 

i-o/j 

0.9/1 

1 o + °- 
1 ■'■'-0.2 

0.9/H 

23 



W,flat,stellar (km S ) 

Fig. 10. Comparison of V Cf i at obtained from stellar and gas kine¬ 
matics for 11 galaxies of our sample 


of the TW method does not result in any late-type galaxy with 
slow bars as predicted by Rautiainen et al. (2008). 

7. Discussion 

7.1. Uncertainties in determiningf the bar pattern speed. 

The pattern-speed measurements using the TW method are sen¬ 
sitive to the determination of the PA of the galaxies. Debattista 
(2002) tested this dependence using N-body numerical simula¬ 
tions. He found that, for PA errors of about 5°, the scatter in 
R is -0.44. The PA of our galaxies was measured by averag¬ 
ing the PA of the outermost isophotes fitted on SDSS r-band 
images. The CALIFA values of i and PA are consistent within 
the errors with those reported in the literature. In particular, the 
median value of the absolute difference between the inclina¬ 
tion of the galaxies reported by LEDA Q dPaturel et al. L120031) 
and our values is A i = \Ileda ~ icalifa\ = 2.6°. There is a 
galaxy (NGC 6945) with a large difference between PA califa 
and PA leda- For this galaxy the value of the PA reported by 
LEDA is 32° ± 66°, which is 94° different from our value. A 
simple visual inspection of the r-band galaxy image showed in 
Fig. [2] ruled out the LEDA PA value. The median value is, ex¬ 
cluding NGC 6945, A PA = \PA leda - PA CALI fa\ = 4.0°. 

Late-type galaxies have strong structures in the outer regions 
of the disc (spiral arms, rings, etc.). These structures can affect 
the orientation of the outer isophotes and subsequently the PA 
value of the galaxy. Thus, it could be that the real PA of the 
discs of our galaxies could be different (more than 3<x) to what 
is obtained from averaging the outermost isophotes. We checked 
this possibility by comparing our photometric PA with the one 
obtained from the symmetrisation of the stellar velocity maps 
(see for more details Falcon-Barroso et al., 2014, in preparation). 
The photometric and kinematic PA of the galaxies agree within 
the 3cr errors for all cases where the velocity map extends to the 
outermost regions of the galaxies where we computed their PAs. 

One of the advantages of integral-field spectroscopy over 
standard long-slit is that the uncertainties in the PA of the galax¬ 
ies can be taken into account by computing Ob along different 
directions according to the errors in PA. This is what we did in 

2 We acknowledge the usage of the HyperLeda database 

(http://leda.univ-lyonl.fr I 


the present work. In particular, we derived Qb in a Monte Carlo 
fashion by varying the PA in the range PA+APA, where APA is 
the uncertainty of the PA given in Table Q] The value of Ob sin/ 
reported in Table [3] corresponds to the mean and dispersion of 
the values obtained for the different PA values. 

Another uncertainty in determining R could be related to the 
determination of V c ji at , which was obtained after a correction 
for asymmetric drift to the observed stellar velocities. This is 
needed due to the stars in the disc of the galaxies do not fol¬ 
low circular velocities. The asymmetric drift correction applied 
in Sect. 6.2 assumes that the galaxies have thin stellar discs. This 
approximation could be unreliable for some of the galaxies. In 
contrast to stars, the gas is almost in circular orbits in galaxy 
discs. We have the gas velocity curve along the major axis for a 
total of 11 of our galaxies (see details in Garcla-Lorenzo et al. 
(2015). FigureflOlshows the comparison between V c ji„, obtained 
from gas and stellar kinematics. For all galaxies except one, 
NGC6941, the agreement between the two quantities is within 
the lcr error bars. This galaxy shows V c ji aUgas IV c ,fiat„ stellar ~ 1-4. 
The values of R for this galaxy would be in the range 0.8 - 1.5 
if we use V c ji al from gas kinematics. 

The convergence of the kinematic integrals could also be an¬ 
other source of uncertainty in the determination of the pattern 
speed. This is especially true for computing the kinematic in¬ 
tegrals using Methods 1 and 3. In these methods we obtained 
these integrals by computing the integrals in Eq. (1) numerically 
by using the binned velocity maps derived from the CALIFA 
datacubes. We estimated the uncertainties of these integrals by 
computing the dispersion of the radial kinematic integrals at 
large radius (beyond the bar radius). For most of our galaxies 
(10/15), these new uncertainties will increase by less than 25% 
the reported uncertainties in Table[4] Exceptions are the galaxies 
NGC5406, NGC6497, NGC6945, UGC03253, and UGC12185. 
For these galaxies the variation on the errors is between 30% and 
35%. 

7.2. Ultralast bars 

Theoretically, the b ar cannot end beyond the corotation radiu s 
of the galaxies (see Contop oulos L Il980t lAt hanassoula , 1992), 
because the main family of orbits constituting the bar (the x\ 
orbits) is unstable beyond CR. This implies that R cannot take 
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Fig. 9. R value as a function of the Hubble type for all galaxies with pattern speed measured by the TW method. The grey points 
represent galaxies from the literature compiled by Corsini (2011), which excludes the sample of Rautiainen et al. (2008), which 
we indicate by the blue rectangles. The red symbols represent the R values of our CALIFA sample. The green stars represent the 
mean values of R for SBO-SBO/a, SBa-SBab, and SBb-SBbc galaxies with O/, determined by the TW method. The blue rectangles 
represent the mean and dispersion (width of the rectangle) values of R from Rautiainen et al. (2008). For clarity reasons, R values 
have been slightly shifted with respect to its Hubble type. 
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values that are less than unity. One of our galaxies (NGC 6497) 
shows R < 1 within the errors (lcr) for our four different mea¬ 
surements. In addition, NGC 0036 and NGC 5205 show 'R val¬ 
ues less than unity within the errors for three of the methods. A 
visual inspection of these galaxies in Fig. [2] reveals no special 
features. There are several ways to explain these R < 1 val¬ 
ues: i) uncertainties associated to the bar and/or corotation radii 
determination; ii) multiple pattern speeds present in the galaxy; 
iii) some other reasons for the TW not being applicable to these 
galaxies. However, we cannot rule out the possibility that they 
might really be ultrafast bars. 

The parameter R is a ratio between the corotation radius and 
the bar radius. Thus, R can be less than one because the bar 
radius is overestimated and/or the corotation radius is underes¬ 
timated. We discarded the case of large bars after a careful vi¬ 
sual inspection of the images. Moreover, these three galaxies 
show values of R below one even if we use the lowest values 
of their bar radius given in Table [2] We computed Rqr using 
the asymptotic flat rotation velocity obtained from the asymmet¬ 
ric drift correction to the LOS velocities of the galaxies. This 
was obtained after some assumptions on the velocity curve of 
the galaxies (see Sect. 5.2). We have found the maximum cir¬ 
cular velocity measured from HI data for one of the ultrafast 


bars (NGC0036 lTheureau et al. Lll998i) . For this galaxy, we have 
Vc.Hi/Vcjiat = 1.14. This means that if assuming V c ,hi the coro¬ 
tation radius and R would be 14% greater than the values given 
in Tables [3] and [4] Thus, if using V c ,hi we would have R\ = 1.0, 
R 2 = 0.7, Ri = 0.7, and R 4 = 0.7. Assuming the same errors as 
given in Tabled this galaxy would only be ultrafast for two mea¬ 
surements. This means that the reason that this galaxy is ultrafast 
could be related to uncertainties associated to the corotation ra¬ 
dius determination. 

We have not found in the literature any HI velocity mea¬ 
surements for the two other ultrafast bars (NGC5205 and 
NGC6497). However, we did recompute the Rqr for NGC5205 
and NGC6497 using the angular velocity curve obtained from 
the gas velocity curve along the majo r axis of these two galaxies 
(see for details iGarcfa-Lorenzo et al7ll2015l) . In these two cases, 
the new values of R were the same as those obtained from the 
stellar kinematics. Thus we cannot discard that these two galax¬ 
ies R might be smaller than one within the errors. 

The observed distribution of R values is given by R 0 b s - 
Rim (X) R err , where 'R„i, s , R,,,,. and R,,,, are the observed, intrin¬ 
sic, and error distributions, respectively. To avoid ultrafast bars, 
we can assume that Rj n , is always larger than 1.0. We can ask 
whether the observed ultrafast bars are only caused by the con- 
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volution of the intrinsic R distribution, where the observational 
errors or additional errors in the TW method are indeed needed. 

The distribution of R\ values in our sample of galaxies peaks 
at R — 1.3 with 37% of the galaxies with R\ < 1. We assume 
that R,„, is centred at 1.3 with dispersion cr ( „, =0.1, in order to 
avoid galaxies with R mt < 1 . The convolution of R mt with the 
observational errors produces a distribution where cr = 0.3 and 
15% of the galaxies show R < 1. Thus, the full fraction of ob¬ 
served ultrafast bars cannot be explained by the convolution of 
R mt with the observational errors. We have to convolve R„„ with 
a distribution R err with cr = 0.5 in order to obtain the 37% of 
ultrafast bars. The dispersion of the new error distribution was 
obtained by adding the dispersion of the observational errors in 
quadrature plus the dispersion of another distribution of addi¬ 
tional errors in the TW method. The additional errors in the TW 
turned out to have a dispersion that is 1.3 greater than the dis¬ 
persion of the distribution of the observed ones. The large extra 
errors that are needed make it unlikely that ultrafast bars are just 
a consequence of new errors not being taken into account in the 
TW method. Similar results were obtained for R 2 , R 3 , and R\ 
distributions. 


The presence of dust lanes within the bar region can affect 
the value of the pattern _speed of bars determined by the TW 
method (see iGerssen & Debattista L l2007t) . This variation de¬ 
pends on the angle that the slits of the TW method cut to the 
bar (APAbar). Thus, dust lanes on the leading edges of a bar tend 
to increase the TW-derived value of Q,, when A PA > 0, and de¬ 
crease it when APAb ar < 0. In addition, this variation depends on 
the amount of extinction. Thus, the variation in the TW-derived 
pattern speed of the bar is about 8 % to 25 % for galaxies with 
Ay ~ 3 in the bar dust lines. This variation increases to 20 % 
40 % for Ay » 8. Two of the ultrafast bars, NGC0036 and 
NGC5205, show APA < 0. This means that for these two galax¬ 
ies, the presence of dust lanes within the bar radius should lead 
to lower TW values of Sib than would the real one. Therefore the 
TW-derived value of R would be greater than the real ones ow¬ 
ing to the presence of dust lanes. We can conclude that the low 
values of R obtained for these galaxies are not due to the strong 
dust lanes within the bar regions. In contrast, A PAb ar > 0 for 
NGC6497. Thus, according to Gerssen & Debattista (2007), the 
TW method will show higher values of Sib and lower values of R 
than the real ones owing to dust lanes within the bar region. For 
this galaxy we obtained the extinction map from the fit of its stel- 
lar population (see lPerez e t al. L 20131: iGonzalez* D elgado et al. L 


popu_ _ ■ _ 

2014; Cid Fernand es et al. . 20131) . The mean and maximum Ay 
values within the bar region of NGC6947 are 0.16 and 0.66, re¬ 
spectively. Thus, we expect maximum uncertainties of 25% in 
the value of R (see Debattista & Gerssen 2007). These uncer¬ 
tainties are not enough to give R > 1 within the errors (see Table 

4). 


Ultra-fast bars have also been observed in other studies of 
barred galaxies. In particular, some of the barred galaxies in the 
sample by Buta & Zhang (2009) show R < 1. These authors 
argue that some of them could be true ultrafast bars and not arte¬ 
facts due to wrong measurements. Other TW measurements of 
bar pattern speed have also shown R < 1 (see iCorsini et al. I 
120071) . The results for the mean value of R and the dependence 
of R on Hubble type presented here do not change if we exclude 
these two galaxies. 

We cannot discard that the TW method cannot be applied 
for the two ultrafast bars owing to these galaxies not following 
some of the assumptions of the method. Further research on this 
direction is therefore needed. 


7.3. Astrophysical implications of fast bars 


Numerical simulations have shown that the bar-pattern speed 
depends on the bar formation and evolution. If a bar forms 
by a global bar in stab ility, then it tends to be a fast 
bar (e.g., ISellwoodl 1198 ll) . Nevertheless, some simulations 
have _shown that bars formed by a gradual bar growth 
( Lvnden-Bell L Il979l). b v interactions with other galaxies 
(Miwa & Noguchi 1 19981), or when the initial bulge-to-disc ra- 
tio is low dCombes & Elmegreen , 119931) can form as slow bars. 
In particular. Combes & Elmegreen (1993) show that late-type 
galaxies have a low-mass concentration at G - k/2 (with k the 
epiciclic frequency) so that at the beginning of bar formation the 
CR is far out in the disc, beyond the disc scale length. In con¬ 
trast, the bar length is determined by the disc scale length. Then, 
the bar starts out as slow. This does not happen for early-type 
galaxies where the CR is shorter and determines the bar length, 
so that the bar is fast from the beginning. According to these 
simulations, and assuming no change in the pattern speed after 
bar formation, we should expect that the bar pattern speeds of 
early- and late-type barred galaxies should be different, because 
they are slower for late-type ones. 

Previous observational works show contradictory results. 
Aguerri et al. (1998) used a sample of ten barred galaxies to 
find a hint of variation in the bar-pattern speed with the mor¬ 
phological type. In the same direction, Rautiainen et al. (2008) 
conclude that late-type galaxies have slow bars, while early-type 
galaxies do not. In contrast, Buta & Zhang (2009) find that there 
is no dependence of R on the morphological type. Before the 
present work, only the pattern speed of three late-type galax- 
ies was obtained by the TW method (e.g.. IGerssen et alTL 12003b 
iTreuth ardt et al. , 2007). These three galaxies turned out to be 
compatible with fast rotators. Our result confirms the previous 
findings, and we clearly show no trend of R with the Hubble 
type. Bars located in early- and late-type galaxies are compati¬ 
ble with being fast bars. This result conflicts with the different 
origins of early and late-type bars proposed in the Combes & 
Elmegreen (1993) simulations. 

The pattern speed of the bars not only depends on their ini¬ 
tial values, but the subsequent bar evolution also can largely 
change the initial pattern speed. In this sense, several numeri¬ 
cal simulations have shown that the interaction between the bar 
and the dark matter halo of the galaxy produce an exchange of 
angular momentum at the galaxy resonances that can change 
the bar-pattern spee d (see i Debattista & Sellwood L fl998l [2000: 
lAthanassoula et al~L 120031: iMartinez-Valpuesta et al. L 2006). 
Our measurements of fast bars across the Hubble sequence in¬ 
dicate that if bars are initially fast, they interchange few angu¬ 
lar momentum with the dark matter halo independently of the 
Hubble type. Therefore bars are fast during their lifetimes (see 
also lPerez et al. LI20121 ). 

Bar formation and evolution are complex and multiparamet- 
ric problems involving several galaxy parameters, such as gas, 
mass fraction, t riaxiality of dark ma t ter haloes, an d disc velocity 
dispersion (see Athanassoul a et al71 12003L120131) . Our observa¬ 
tional results are important for constraining the different forma¬ 
tion and evolution scenarios proposed by future numerical sim¬ 
ulations. 


8. Conclusions 

We have used integral field spectroscopy data for 15 CALIFA 
barred galaxies to determine their bar-pattern speed using the 
model-independent TW method. Both the FOV and the spatial 
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resolution of the CALIFA data make this dataset appropirate for 
studies of extended disc galaxies, as proposed here. Integral-field 
data have the advange of solving problems such as the centring 
and the disc PA uncertainties that can affect the measured value 
of the bar-pattern speed by the TW method when long-slit spec¬ 
troscopic observations are used. 

The sample of galaxies selected for this study spans a wide 
range of morphological types from SBO to SBbc with most of 
them being spiral galaxies (SBb-SBbc). This sample fills the gap 
present in the literature around spiral galaxies for which the TW 
has not been extensively applied. 

For each galaxy we determined their bar pattern speed us¬ 
ing four different approaches to the TW method. The kinematic 
integrals of the TW were obtained directly from the stellar ve¬ 
locity maps, after summing up the spectra along the correspond¬ 
ing pseudo slits and then measuring the velocity. In addition, 
we used light and mass as weights of the kinematic and photo¬ 
metric TW integrals. We also obtained the distance-independent 
quantity R = /^crAa for each galaxy. The corotation radius of 
the galaxies was approximated by Rcb. = V c ,flat/ ^b. where V Ci fl at 
is the asymptotic circular velocity of the galaxies derived after 
the asymmetric drift correction was applied to the LOS velocity 
along the galaxy’s major axis. 

Our CALIFA sample and those galaxies available in the lit¬ 
erature with bar pattern-speed measurements by the TW method 
make a full sample of 32 galaxies from SBO to SBbc morpholog¬ 
ical types. This full sample is appropriate to studying the possi¬ 
ble dependence of the bar pattern speed on the morphological 
type. 

We found that the mean value of R for all galaxies of the 
sample is in the interval 1.0-1.1. We computed, through Monte 
Carlo simulations, that we cannot rule out (at 95% level) the 
presence of a fast bar (R < 1.4) in any of our galaxies. The mean 
value of R found for our sample agrees with hydrodynamic al nu¬ 
merical simulations of bar galaxies that predict R in the interval 
[1.0,1.4]. We did not find any trend of R with the Hubble type, 
so early and late types are both compatible with being fast bars. 

In the near future, we will increase the number of pattern- 
speed measurements by the TW method using more barred 
galaxies observed in the CALIFA sample. This will be useful for 
enlarging the sample of late-type galaxies, decreasing the uncer¬ 
tainties in the mean values of R, and analysing eventual trends 
in the bar-pattern speed with other galaxy parameters. 
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